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I.  INTRODUCTION 


Th?.  ever  increasing  ability  of  weapons  delivery  systems  to  pinpoint 
their  target  has  prompted  survivability  analysts  and  designers  to  upgrade 
the  resistance  capacity  of  protective  structures.  A  hardened  reinforced 
concrete  structure  covered  by  protective  layer(s)  of  soil,  concrete 
burster  slab,  rock  rubble,  etc.,  is  required  to  provide  protection  against 
the  very  high  overpressure  from  the  near  field,  high  energy  explosion 
of  conventional  (non-nuclear)  weapons.  Conventional  weapons  detonated 
near  a  shallow-buried  structure  result  in  a  rapidly  attenuating,  short 
duration  pressure  pulse,  which  is  in  sharp  contrast  to  the  well  documented 
long  duration  shock  front  associated  with  nuclear  explosions.  The  in¬ 
herent  strength,  stiffness  and  stability  of  boxlike  structures  makes 
them  suitable  for  use  in  such  a  system.  This  project  is  concerned 
with  the  dynamic  response  of  reinforced  concrete  one-way  elements,  the 
majo**  structural  element  of  boxlike  structures. 

Recent  experimental  testing  has  been  performed  in  an  effort 
to  better  define  the  characteristics  of  underground,  non-nuclear  explosions, 
and  their  effects  on  underground  structures.  These  model  tests  provide 
much  needed  information  about  structural  failure  of  reinforced  concrete 
slabs  and  the  in-structure  environment  to  which  the  personnel,  equipment, 
weapons,  etc.  would  be  subjected. 

The  results  of  these  model  tests  of  reinforced  concrete  load 
resisting  elements,  summarized  in  Appendix  A,  indicate  that  the  range 
of  the  failure  modes  vary  from  the  conventional  plastic  hinge  mechanism 
(standard  static  response  engineering  procedure)  to  material  degradation 


due  to  direct  intensity  of  the  blast  pressure  (no  structural  member 
response  involved).  The  present  work  is  primarily  directed  toward  the 
investigation  of  an  approximate  analytical  model  cf  the  failure  modes 
associated  with  structural  member  response.  Available  techniques  for 
theoretically  determining  direct  material  degradation  characteristics 
(i.e. ,  spalling)  will  also  be  discussed. 

Analytical  and  experimental  work  in  the  aeronautical  field  for 
years  has  been  devoted  to  the  response  of  metal  beams  and  plates  with 
solid  rectangular  cross-sections.  Dynamic  rigid-plastic  theory  has 
been  developed  for  the  high  intensity  loading  associated  with  the  blast 
phenomena.  This  approach  is  reasonable  when  the  energy  input  to  the 
system  is  significantly  larger  than  the  maximum  "elastic"  strain  energy; 
in  other  words,  the  elastic  deformations  are  negligible  compared  to  the 
plastic  deformations. 

Simply  supported  and  fixed  ended  beams  without  axial  restraint, 

subjected  to  uniformly  distributed  impulsive  loading,  have  provided  the 

basis  for  many  investigations.  Recently,  the  temporal  and  spatial 

limitations  for  the  loading  required  by  these  closed  form  solutions 
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was  overcome  by  the  development  of  a  computer  program,  CONCRE  ’  , 
which  computes  the  deflection  of  the  center  of  a  reinforced  concrete 
beam  or  plate  caused  by  a  linear  or  blast  load.  Reference  5  provides 
experimental  evidence  that  buried  fixed  ended  reinforced  beams  with  no 
axial  restraint  do  exhibit  a  conventional  plastic  hinge  failure  mode, 
demonstrating  that  under-reinforced  concrete  elements  possess  the 
necessary  ductility  for  its  response  to  be  idealized  by  a  rigid,  perfectly 


plastic  constitutive  relationship.  If  the  blast  pressure  is  sufficient, 
the  compressive  reinforcement,  in  addition  to  the  tensile  reinforcement, 
may  fracture  along  well  defined  yield  lines,  resulting  in  total 
disengagement.  The  sections  between  hinges  essentially  remained  intact. 
Because  of  the  lack  of  in-plane  translation  restraint  on  these  one-way 
slabs,  the  potential  membrane  capacity  of  the  member  is  not  utilized 
as  it  might  be  in  a  real  boxlike  structure. 

When  the  edges,  or  at  least  the  corners,  of  a  transversely  loaded 
member  are  restrained  against  in-plane  motion,  lateral  deflections 
beyond  a  certain  level  (a  >_  0.3t)  are  accompanied  by  stretching  of 
the  middle  surface.^  As  the  magnitude  of  the  maximum  deflection  in¬ 
creases,  the  resulting  membrane  force  can  even  predominate  in  providing 
the  lateral  load  resistance  of  an  axially  restrained  member.  Test 
results11  Indicate  that  rectangular  cross-section  steel  beams  with 
full  end  fixity  (flexural  plus  axial)  may  develop  a  static  load  capacity 
six  or  seven  times  the  load  carrying  capacity  of  a  freely  supported 

beam.  Visual  examination  of  the  photographic  records  of  laboratory  scale 

12 

blast  loading  tests  by  Abrahamson  on  rectangular  cross-section  metal 
beams  and  plates  show  maximum  deflections  that  are  multiples  of  the 
member  thickness.  For  the  experimental  results  of  blast  loads  on 
reinforced  concrete  models  previously  referenced,  distinctive  yield  lines 
(plastic  hinges)  are  not  generally  observed  to  form  in  the  test  specimens 
which  are  provided  in-plane  end  restraint  against  translation  by  their 
supports. 

Much  of  the  more  recent  analytical  efforts  have  been  devoted  to 
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describing  the  response  of  reinforced  concrete  elements  considering 
only  the  flexural,  exclusive  of  the  extension,  strain  energy  in  their 
formulation.  The  inclusion  of  the  effects  of  an  axial  restraint  in  the 
analysis  not  only  reoiires  additional  terms  in  the  differential  equation 
describing  the  structural  response,  but  also  M,  the  plastic  moment 
capacity  of  the  cross-section,  is  reduced  below  MQ,  the  "full"  plastic 
moment  capacity  of  the  cross-section.  This  capacity  reduction  depends 
on  the  magnitude  of  N,  the  tensile  force  acting  on  the  cross-section. 

If  the  intensity  of  the  loading  is  sufficient  to  cause  N  to  equal  NQ, 
the  "full"  plastic  tensile  capacity  of  the  cross-section,  then  the 
flexural  capacity  M  of  the  cross-section  equals  zero  for  the  combined 
loading.  With  the  assumption  that  the  normal  tensile  force  at  the 
location  of  plastic  flow  is  constant  along  the  member,  the  succeeding 
response  would  resemble  that  of  a  plastic  string  or  cable. 

This  report  summarizes  a  rigid-plastic  dynamic  analysis  procedure 
for  an  axially  unrestrained  beam,  then  transitions  to  an  investigation 
of  an  axially  restrained  member,  which  develop  internal  bending  moment 
and  axial  force.  These  procedures,  previously  developed  for  ductile 
metal  beams,  will  be  applied  to  reinforced  concrete  members. 


Loading  Cases 

Blast  loading  is  characterized  by  a  pulse  with  an  instantaneous 
rise  to  a  peak  positive  pressure  -Followed  by  a  decay  to  zero.  This 
positive  pressure  duration  is  followed  by  a  period  of  negative  pressure 
(suction),  as  illustrated  by  the  insert  in  Figure  1.  The  culmination 
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of  considerable  effort  to  document  blast  load  parameters  is  summarized 

13 

graphically  in  the  well  accepted  reference  AFM88-22  .  An  abbreviated 
copy  of  this  graph  is  provided  in  Figure  1,  where  the  various  blast 
loading  parameters  are  plotted  against  the  scaled  distance  R/W1^3. 

These  values  are  for  free  field  air  blast,  with  modifications  required 
to  account  for  the  reflection  effects  when  encountered  Ly  a  structure. 

A  similar  widely  accepted  procedure  for  underground  blast  parameters 
for  near  field  conventional  explosions  Is  not  available.  At  this  point 
in  the  report  a  brief  description  of  a  recently  publicized  procedure 

for  determining  underground  blast  parameter  values  is  provided. 
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In  recent  publications  by  Southwest  Research  Institute  *  ,  an 

empirical  approach  was  used  to  determine  underground  blast  load  parameters. 
A  more  sophisticated  analytical  approach  was  not  used  because,  among 
other1  reasons,  no  recognized  equation-of-state  are  available  for  different 
soils  exposed  to  ground  shock  from  nearby  detonations.  SRI  used  the  pi 
theorem  to  determine  dimensionless  ratios  (pi  terms)  containing  the 
appropriate  parameters.  Since  only  dimensionless  ratios  have  to  be 
the  same  in  order  for  two  systems  to  be  equivalent,  properly  scaled  model 
tests  were  then  used  to  provide  the  necessary  information  to  determine 
the  functional  format  for  soil  particle  velocity,  U,  and  soil  displacement, 
X,  in  the  free  field.  For  example,  extensive  experimental  data  illus¬ 
trated  in  Figure  2,  was  used  to  determine  the  following  function  for 
the  soil  particle  velocity: 


-4i  -t  -• 
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0.00617  [—5^ 

W ¥> 
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3  mass  density  of  the  soil  or  rock  (lb-sec  /ft  ) 

=  seismic  P-wave  velocity  in  the  soil  or  rock  (ft/sec) 

o 

*  atmospheric  pressure  (Ib/ft  ) 


The  minimum  value  of  scaled  distance  for  this  experimental  data  is 
approximately  2.4,  which  is  sufficient  to  cause  structural  mode  failure 
in  burled  reinforced  concrete  members.  It  is  apparent  that  additional 
experimental  data  is  needed  to  extend  this  function  to  a  lower  limit 
for  scaled  distance. 

Using  the  free  field  soil  particle  velocity  and  the  soil  displace¬ 
ment,  the  free  field  pressure  can  be  approximated.  As  with  obstructions 

in  air,  buried  structures  are  subject  to  the  effects  of  a  reflected 

1 5 

pressure  (amplified  free  field  pressure).  Experimental  evidence  is 
presently  being  generated  in  order  to  evaluate  the  effects  of  the  load- 
deformation  properties  of  a  buried  structural  element  on  the  magnitude 
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of  the  reflected  pressures  felt  by  the  structure.  The  future  utilization 
of  this  research  will  require  a  more  realistic  model  of  the  deformational 
response  of  the  structural  elements  that  comprise  buried  box  structures. 
Presently,  widely  used  procedures  simply  double  the  f^ee  field  pressures 
in  order  to  determine  the  design  loads  on  underground  structures.  This 
implicitly  assumes  that  the  "rigid"  walls  terminate  the  soil  particle 
motion  normal  to  the  plane  of  the  wall. 

Using  a  reflected  pressure  coefficient  of  2.0,  SwRI  published  the 
following  equation  in  reference  14  for  the  determination  of  the  impulse 
applied  at  any  point  x  on  a  buried  wall. 


Similar  equations  are  presented  for  the  peak  reflected  pressure  dis¬ 
tribution.  Except  for  explosions  close  enough  to  result  in  direct  material 
degradation,  a  valid  assumption  that  greatly  simplifies  structural 
response  investigations  considers  blase  loadings  impulsive.  Combined 
with  mass  distribution  information,  the  velocity  imparted  to  a  structura1 
element  can  be  computed.  It  should  be  noted  that  reference  14  suggests 
that  the  mass  of  the  soil  between  the  explosive  charge  and  the  structural 
element  be  included  with  the  mass  of  the  element  in  order  to  obtain  a 
more  accurate  estimate  of  the  overall  dynamic  response  of  the  structural 
element. 


While  no  success  was  had  in  locating  a  published  functional  re¬ 
lationship  for  underground  blast  parameters  as  the  charge  location 
approaches  that  of  a  contact  charge,  this  range  of  scaled  distances 
was  addressed  for  air  blasts  in  1  recent  publication.  Kot  et  al  *s^ 
procedure,  which  incorporates  the  previously  referenced  air  blast  para¬ 
meter  values  from  AFM88-22,  is  applied  to  scale  distances  as  small  as 
0.2,  which  is  well  within  the  deflagration  zone  of  a  spherical  charge. 

The  spall  analysis  for  this  level  of  intensity  of  loading  is  described 
in  a  latter  section  of  this  report.  For  now,  suffice  it  to  point  out 
that  the  procedure  involves  using  direction  cosines  for  the  components 
normal  to  the  member  along  their  inclined  scaled  distance.  The  temporal 
variability  of  these  loadings  has  been  ignored  as  the  analysis  procedure 
assumes  maximum  impulsive  loadings  irrespective  of  the  time  of  arrival 
at  the  various  locations  along  the  member.  Short  of  an  extensive  finite 
element/finite  difference  analysis,  considerable  determination  of  the 
structural  response  of  reinforced  concrete  members  can  be  had  using 
uniformly  distributed  loadings  in  our  analysis. 

In  order  to  obtain  closed  form  solutions  to  blast  response  problems, 
a  considerable  amount  of  the  literature  is  limited  to  uniformly  distributed 
impulsive  transverse  loads  on  flexural  members.  Considering  ■'.he  arrival 
time  required  for  near  field  blast  wave  to  reflect  at  different  locations 
along  the  length  of  a  structural  member,  the  complicated  loading 
distribution  is  certainly  not  conducive  to  representation  by  a  single 
decaying  time  function.  Even  Ross1  computer  program,  CONCRE  ’  ,  which 
is  not  constrained  to  impulsive  loadings,  does  not  permit  variable 
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temporal  functions  at  different  locations  along  a  beam.  It  should  be 
pointed  out  that  spatial  variations  allowed  in  the  procedure  can  result 
in  a  theoretical  response  that  violate  the  inherent  requirement  in 
rigid-plastic  response  that  internal  moments  at  all  locations  along  the 
member  not  exceed  the  plastic  moment  capacity  at  plastic  hinge  locations. 

The  spatial  distribution  of  the  blast  loading  for  the  determination 
of  the  structural  response  of  transversely  loaded  members  in  this  report 
will  be  limited  to  one  that  is  uniformly  distributed.  The  significance 
of  including  membrane  effects  can  be  demonstrated  without  a  more  exact  . 
temporal  and  spatial  definition. 

Member  Strength 

The  reinforcement  of  concrete  with  ductile  steel  has  long  been 
recognized  as  a  means  of  overcoming  the  inherent  tensile  weakness  of 
concrete  when  used  in  a  composite  structural  member.  Whitney's  Method 
for  computing  the  ultimate  flexural  capacity  of  a  cross-section  ignores 
the  tensile  strength  of  the  concrete  and  replaces  the  nonlinear  compressive 
stress  distribution  (see  Figure  3)  with  a  mathematically  equivalent  rec¬ 
tangular  stress  block.  The  same  assumption  is  used  to  generate  interaction 
diagrams  for  use  with  beam-column  member  capacity  determination.  As 
the  name  of  the  member  implies,  a  beam-column  member  is  subjected  to 
compressive  forces  in  addition  to  flexure.  An  illustration  of  a  typical 
interaction  diagram  for  a  member  with  longitudinal  reinforcement  in  two 
parallel  faces  is  provided  in  Figure  3.^  Note  that  Mq  is  the  pure 
flexural  capacity  of  the  member.  In  design  practice,  the  steel  ratio  in 


pure  flexural  equations  refers  only  to  the  quantity  of  tensile  steel, 
while  in  beam-column  design  aids  the  steel  ratio  is  based  on  the  total 
amount  of  longitudinal  steel  in  the  cross-section. 

As  this  work  is  concerned  with  the  relatively  uncommon  influence 

of  an  axial  tension  on  the  response  of  a  flexural  member,  the  typical 

interaction  diagram  in  Figure  3  is  modified  to  include  the  combined 

effects  of  tension  and  bending  on  the  member  capacity  by  extending  the 

ordinate  on  the  negative  (or  tensile)  side  of  the  axis.  NQ,  the  maximum 

value  for  the  tensile  force  when  no  bending  moment  is  present,  is  given 

by  Nq  =  Asfjy*  with  Fy  equal  to  the  yield  strength  of  the  reinforcing  steel. 

Because  the  concrete  will  crack  before  the  steel  yields,  only  the  steel 

is  assumed  to  resist  the  axial  tensile  force  applied  to  the  member  at 

ultimate.  The  neutral  axis  position  for  combined  bending  and  tension 

is  even  smaller  than  its  value  for  pure  bending,  and  when  the  axial  tension 

is  sufficiently  large  in  magnitude,  the  neutral  axis  falls  theoretically 

beyond  the  edge  of  the  cross-section  (there  is  no  compression  stress 

in  the  member).  The  dashed  line  in  Figure  3  connecting  the  limiting 

values  for  axial  tension  NQ  and  bending  moment  MQ  illustrates  an  assumed 

shape  of  the  interaction  diagram  for  this  region.  This  curve  will  be 

approximated  in  this  study  by  the  relationship  (M/MQ)  +  (N/NQ)2  -  1, 

18 

which  is  used  by  Symonds  and  Mentel  to  define  the  relationship  between 
axial  tension  and  bending  moment  for  a  solid  rectangular  ductile  cross- 
section. 

The  modes  of  failure  of  beams  and  one-way  slabs  exhibited  in 
field  tests  with  conventional  blast  loads  range  from  simple  plastic 


hinge  mechanisms  through  direct  material  degradation  due  to  the  intensity 
of  the  compressive  blast  pressure  on  the  exposed  surface  of  the  member. 

The  reason  for  including  the  membrane  behavior  in  the  analysis  is  an 
attempt  to  explain  the  severe  cracking  that  accompanies  the  structural 
response  failure  modes  even  when  the  intensity  of  the  blast  pressure  is 
not  sufficient  to  have  caused  the  damage  by  either  spalling  associated 
with  the  compressive  wave  traveling  through, the  member  thickness  being 
reflected  as  a  tension  stress  at  the  rear  free  surface,  or  direct  com¬ 
pression  blast  pressure  sufficient  to  fracture  the  concrete  in  the 
transverse  direction  due  to  Poisson's  effects.  Also  when  considering 
the  structural  response  of  a  member  following  this  initial  level  of 
direct  material  damage  to  the  concrete,  only  the  membrane  capacity  01 
the  reinforcing  steel  is  available  to  the  member.  The  behavior  of 
concrete  as  a  material  is  particularly  relevant  to  comprehending  the 
different  failure  modes  exhibited  by  the  experimental  tests.  Therefore, 
even  though  concrete's  tensile  strength  is  not  explicitly  considered  for 
design  purposes,  an  understanding  of  concrete  tensile  strength  is  -equired 
to  explain  various  failure  modes. 
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Flexural  Response 


References  18  and  19  provide  concise  documentation  of  the  rigid- 

plastic  response  of  both  simply  supported  beams  and  fixed  ended  beams, 

with  no  axial  constraint,  to  a  uniform  oynamic  load  q(x,t)  =  q(t). 

Consistent  with  these  procedures,  the  effect  of  the  change  of  normal  and 

shear  forces  on  the  ultimate  bending  moment  are  not  considered.  When 

the  effects  of  axial  restraint  are  ignored,  the  rigid  plastic  solution 

is  quite  straightforward.  The  moment-curvature  relationship,  ignoring  elas 

tic  response,  results  in  the  basic  requirement  that  no  rotation  occurs 

if  the  internal  moment  is  less  than  MQ,  the  full  plastic  moment  capacity 

of  the  cross-section.  However,  unlimited  plastic  flow  can  occur  if  the 

upper  limit  moment  is  maintained  at  a  cross-section. 

To  review  briefly,  for  a  uniformly  loaded,  simply  supported  beam, 

the  static  collapse  mechanism  develops  a  plastic  hinge  at  the  beam 

center.  Figure  4a  illustrates  that  the  sum  of  moments  about  the  support 

2 

yields  the  static  collapse  load  as  qg  =  8Mq/L  .  Had  the  load  been 

applied  dynamically,  the  response  would  depend  upon  the  intensity  of 

the  blast  pressure.  Consistent  with  the  finding  of  other  investigators, 

if  the  maximum  intensity  of  the  blast  loading  is  qs  i  q  1  3q$;  then  the 

failure  mechanism  is  identified  as  having  a  single  stationary  hinge  at 

the  center  of  the  beam  span,  x  =  L/2.  Refer  to  mechanism  I  in  Figure  4a. 

A  consideration  of  the  freebody  of  the  motion  results  in  the  following 

19 

equation  for  the  bending  moment  on  the  cross-section  at  x  =  L/2; 

M(Xit) .  Mo(^ .  4 (y  Sit!  ♦  «  -  4*) 

At  Kr 
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and  for  the  equation  of  motion 


8 


of  mechanism  I: 


0  <_  t  <_  T 

t  >  t 


For  the  mechanism  to  be  valid,  both  halves  of  the  beam  must  behave  as 
rigid  plates,  M(x,t)  <  MQ. 

For  more  intense  dynamic  loading,  q  >  3q  ,  the  structural  response 
deformation  pattern  in  the  first  phase  of  motion  is  given  in  Figure  4b. 
Note  the  two  symmetrically  located  non-stat ionary  plastic  hinges.  The 
center  portion  between  the  plastic  hinges  translates  vertically  with  a 
linear  velocity,  while  the  two  rigid  end  segments  rotate- about  their 
supports.  Expressing  the  continuity  of  velocity  at  the  intersection 
of  the  end  segment  with  the  horizontal  central  section  results  in  the 
following: 


wa(t)  =  V 

Mathematically,  the  boundary  between  the  segments  changes  its  position 

during  intense  dynamic  loading,  resulting  in  a  bending  wave  propagation 

horizontally  along  the  beam  (or  a  moving  yield  line).  The  following 

19 

equation  from  Henrych's  work  shows  that  the  boundaries  are  moving  in 
a  direction  from  the  supports  towards  the  center  of  the  beam: 


EJTt 

a(t)  =\Hr - 

\J/Jq(t)dt 


At  time  t  =  t-j,  a(t-j)  =  L/2;  therefore,  the 
a  single  hinge  at  the  center  of  the  member. 


two  hinges  are  united  into 
The  second  phase  motion. 
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t  >  t-j,  is  then  simply  that  of  a  stationary  plastic  hinge  at  the  center 
of  the  beam,  as  shown  previously  in  Figure  4a.  The  time  of  maximum 
deflection,  tm»  is  calculated  from  the  boundary  condition  of  zero 
velocity.  Therefore,  for  a  blast  intensity  sufficient  to  initiate  a 
mechanism  II  response,  the  value  for  the  maximum  deflection  at  midspan 
of  the  member  follows  directly  as  the  addition  of  the  first  and  second 
phase  response:  vmax  =  v-j(t-j)  +  Vg’(t  "  t-j).  The  computer  program, 
CONCRE,  utilizing  numerical  integration  facilitates  the  process  of 
actually  performing  the  calculations  involved  with  the  various  pressure  ' 
decaying  functions  that  are  common  to  blast  loading  problems.  This 
computer  code  also  permits  spatial  loading  variations;  however,  moments 
exceeding  the  plastic  capacity  of  the  cross-section  MQ  may  occur  between 
hinges.  In  Appendix  B  these  equations  of  motion  are  derived  using  the 
Lagrange  equation. 

For  loadings  that  can  be  considered  impulsive,  /£q(t)dt  =  i, 
considerable  simplification  can  be  achieved.  Consult  reference  19  for 

the  closed  form  solutions  of  impulsively  loaded  beams,  i.e., 

P  _  1  m£^V^ 
rinal  '  3  ~R 

o 

Beams  with  Axial  Constraints 

18 

Symonds  and  MenteVs  rigid-plastic  solution  procedure  for  the 
response  of  a  beam  with  axial  constraints  assumes  that  the  normal  force  N 
on  a  cross-section  is  constant  along  the  member  (see  Figure  5).  Plastic 
deformations  are  assumed  to  occur  when  M  and  N  satisfy  the  plasticity 
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criteria  that  (M/Mq)  +  (N/NQ)2  =  1. 

The  flow  rule  is  used  to  relate  the  plastic  deformations  that  result 
when  this  criteria  is  satisfied  recognizes  that  the  strain  rate  vector 
is  normal  to  the  yield  curve,  which  is  expressed  mathematically  as 

V  =  -  N 

MQ$  No 

where  e  is  the  extension  strain  and  ^  is  the  curvature,  and  the  dot 
indicates  time  derivatives. 

Considering  an  intensity  of  loading  sufficient  to  cause  an  initial 
plastic  hinge  location  not  at  the  center  of  the  member,  four  unknowns, 

M,  N,  a,  and  u>,  will  need  to  be  determined.  The  previously  described 
equations  for  continuity  of  velocity  at  hinge  locations  and  the  flow 
rule  relating  M  and  N  provide  two  of  the  four  equations  needed.  The 
addition  of  the  axial  force  effect  on  the  previously  examined  angular 
acceleration  equation  yields 

■j  ma3u)  =  -  M  -  N  6 

where  6  is  the  deflection  of  the  center  section  of  the  member.  An 
additional  condition  resulting  from  the  axial  constraint  is  formulated 
from  a  consideration  of  the  rotation  of  segment  AB  as  it  moves  through 
an  angle  A0  in  a  small  time  interval  At.  This  response  causes  the 
plastic  hinge  to  move  to  $  with  the  length  of  the  segment  increase 
by  ax.  The  curvative  change  Tp  in  the  segment  Ax  is  given  by  Ae/Ax, 
and  the  geometrical  requirement  in  Figure  5  is  eAx  =  (Ae)<5. 

As  formulated  by  Symonds  and  Mentel ,  this  relation  between  extensional 
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strain  and  curvature  change,  combined  with  the  flow  rule  provide  the 
following  equation: 

N_  V 
n-m; 

During  the  time  the  hinges  are  moving  toward  the  center  of  the 
member;  <S  =  Vt.  Therefore,  the  initial  phase  response  of  mechanism  II 
behavior  is  found  from  the  solution  of  simultaneous  equations.  The 
result  of  this  is: 

mV3 

for  which  the  solution  of  this  differential  equation  is: 


/  ,/  N0w  v 


This  initial  phase  motion  terminates  at  time  tj,  when  a(t)  =  with 


8  V,  Therefore 


+  _  mV* 
‘1-6H7 


'V  V 

st/ 


The  two  halves  of  the  beams  (rigid  bars)  then  rotate  about  the  supports 
during  the  second  phase.  The  permanent  curve  change  and  axial  extensions 
are  assumed  to  occur  at  the  plastic  hinge  at  midlength  of  the  beam. 

The  unkiiown  quantities  are  M,  N,  w  =  e,  and  6.  The  same  four  simultaneous 
equations  as  for  the  initial  phase  response  are  solved  with  a(t)  =  r 
(plastic  hinge  at  midspan).  The  solution  of  the  resulting  differential 
equation,  and  simplification  of  the  results  in  reference  18  furnishes 
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wal 


the  following  equation  for  the  final  deflection  6. 


N  P 

where  the  dimensionless  parameter  e  =  o 

1 2M_ 


and  a  =  m  z\l 

6fT 

o 

For  beams  with  fixed  ends,  the  final  deflections  can  be  determined 
simply  by  replacing  MQ  with  2MQ  in  the  equations  for  final  deflections. 
This  rather  simple  analysis  provides  results  for  the  effect  of  axial 
restraint  on  the  transverse  response  of  beams  assuming  rigid  plastic 
behavior.  However,  the  foregoing  procedure  is  valid  only  so  long  as 
N/Nq  <1.  At  succeeding  intervals  of  time,  the  axial  force  probably 
remains  constant  and  equal  to  NQ,  which  establishes  the  bending  moment 
is  zero  throughout  the  beam.  The  classical  partial  differential 
equation  describing  the  future  behavior  of  the  member  as  a  plastic 


string  is: 


2  2 
a_y  _  m  ry 


For  the  details  of  the  derivation  of  the  response  of  a  member  in 
the  plastic  string  phase  you  are  referred  to  section  4  of  reference  18. 
The  approximate  conservative  solution  from  this  publication  is  given  by 


the  following  equation 


si  < 

Z  1  N  M 
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The  following  example  problem  is  used  to  illustrate  the  relative 
magnitude  of  deflection  when  axial  constraint  is  considered  in  the  analysis 


Simply  Supported  Concrete  Beam  with 
Symmetrical  Tensile  and  Compressive  Steel 

Span  =  L  =  80  in.  fc  =  6.0  ksi 

Depth  =  t  =  6  in.  f  =  60.0  ksi 

Width  =  b  =  12  in. 

Steel  Ratio  =  pt  =  0.01 

Cover  for 

Steel  =  0.75  in. 

Ratio  of  distance  between  centroid  of  outer  rows  of 

bars  and  the  depth  of  the  cross-section  -  y=  0.75 

Properties  of  cross-section: 

m  *  mass  per  unit  length 

.  (12  In.  x  6  in. HO. 083  Win3)  „  1  kip 
386  in/sec2  1000  1b- 

2 

m  =  15.5  x  lO"6  -- r-s—  per  inch  length  of  member 
Nq  =  Axial  Yield  Capacity 

=  Asfy  =  ptAgfy  =  (0.01  x  72  in2) (60  K/in2) 

Nq  =  43.2  kips 

Mq  =  Pure  Moment  Capacity  of  Cross-section 
Using  the  ACI  interaction  diagrams^,  the  pure  flexural  capacity 
for  this  doubly  reinforced  cross-section  is  determined  from  the  following 
equation  and  appropriate  moment  coefficient: 


i  Mo 

V 


0.2331 
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For  a  capacity  reduction  factor  0  =  1,  the  flexural  capacity 
MQ  =  100.7  in-k. 

For  an  initial  velocity  VQ  equal  to  500  in/sec  imparted  to  the 
beam  by  an  impulsive  blast  load, 
with  no  axial  constraint: 

6-  mil!/2 

T  _  0 

l  “  3M_ 

0  2 

.  (15.5  x  10“6  ^=SSC  )(40  in)(500  in/sec)2 


*  0.51 


with  axial  constraint: 


*  °|n0 


500  in/sec 


15.5  x  10 


-6  k-sec‘ 


100.7  in-k 


Of 

4  .  0.24 

Figure  6  illustrates  the  effects  of  the  initial  velocity  on  the 
relative  magnitude  of  the  final  midspan  deflection  when  axial  restraint 
is  considered.  The  curve  demonstrates  that  for  a  typical  reinforced 
concrete  member  the  deflections  calculated  for  a  simply  supported 
unrestrained  beam  will  be  much  larger  than  the  deflections  determined 
for  an  axially  restrained  beam.  Even  though  these  results  are  based  on  an 
approximate  procedure,  it  is  certainly  appropriate  for  design  ptucedures. 
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In  order  to  utilize  the  previous  information,  a  failure  criteria 
must  be  determined  to  relate  the  rupture  of  the  reinforcing  steel  to  the 


midspan  deflection  <sf.  Fran  standard  assumptions  relative  to  membranes, 
the  following  equation  for  the  longitudinal  strain  in  the  reinforcing 


steel  is: 


Rearranging:  <5^  2 

%  IT 

From  ASTM  A615,  the  tensile  specifications  for  Grade  40  and  Grade  60 
reinforcing  steel  require  material  elongations  (ex)max  'from  8%  to 
for  the  size  reinforcing  bars  normally  used  in  blast  resistant  structures. 
For  this  example  problem  the  steel  would  theoretically  rupture  in  tt;e 
axially  restrained  member  at  {5^/Jt)  values  of  0.255  and  0.298,  respectively. 


Spalling  of  Concrete  Walls 

When  computing  the  structural  response  of  reinforced  concrete 
elements  to  intense  blast,  authors  have  invariably  ignored  the  loss  of 
the  composite  properties  of  the  cross-section  due  to  spallation,  etc. 

While  the  structural  response  is  generally  considered  as  a  reaction 
to  impulsive  loads,  several  internal  stress  wave  reflections  associated 
with  spalling  have  occurred  before  the  structural  motion  begins.  As 
a  consequence,  the  small  sized,  high  velocity  spall  debris  (normally 
considered  to  be  comprised  of  the  concrete  cover  over  the  reinforcement, 
on  the  back  side  of  the  wall)  has  been  discharged  before  the  overall 
structural  motion  has  commenced.  As  stated  in  reference  16,  the  smaller 
velocity  imparted  to  the  larger  spall  particles  causes  them  to  stay  with 
the  structural  element  and  to  be  ejected  later  by  the  larger  wall 
velocities  associated  with  structural  response. 

The  immediately  occurring  extensive  internal  damage  to  the  concrete 
associated  with  the  internal  tension  stresses  resulting  from  the  reflected 
compressive  blast  waves  would  most  certainly  eliminate  the  flexural  cap¬ 
acity  of  the  member.  The  structural  motion  resistance  would  then  be 
limited  to  the  membrane  capacity  of  the  reinforcing  steel. 

Literature  documenting  an  analytical  procedure  for  determining 
spallation  of  underground  walls  was  not  found;  however,  comprehensive 
references  on  spallation  as  the  result  of  air  blast  on  a  structural  wall 
were  discovered.  The  principles  applied  to  these  walls  apply  equally 


well  to  a  buried  wall,  with  the  exception  that  the  relationship  between 
charge  weight,  standoff  distance,  and  the  resulting  load  parameters 


( i . e. ,  pressure,  impulse,  etc.)  would  be  different.  For  both  situations, 
identical  loading  parameters  on  a  wall  would  yield  equivalent  responses. 

The  analytical  work  done  on  this  subject  that  is  referenced  here 
required  considerable  simplification  of  a  very  complex  problem.  The 
basis  for  the  simplified  analysis  is  well  known  and  is  simply  restated 
here  from  reference  16.  “The  phenomena  of  spall  in  brittle  materials 
occurs  when  strong  tension  waves  are  reflected  into  the  wall  from  the 
free  surface  at  the  back  face  and  interact  with  the  decaying  compression 
wave  in  such  a  way  as  to  produce  locally  tension  stresses  which  exceed  . 
the  dynamic  tensile  rupture  strength  of  the  material."  An  estimate 
of  the  threshold  was  proposed  in  reference  14  using  the  following 
assumptions:  (1)  uniform  loading  of  wall,  (2)  triangular  pressure¬ 
time  history,  which  directly  relates  the  peak  reflected  pressure  Pr> 
reflected  positive  specific  impulse  ir,  and  positive  pulse  duration  T, 

(3)  no  attenuation  of  the  pressure  wave  through  the  wall,  and  (4)  a 
stress  wave  velocity  through  the  wail  given  by  v  =  &  ,  where  E  is 
Young's  Modulus  and  p  is  density.  The  results  of  this  analysis  for  pre¬ 
dicting  the  threshold  of  spallation  from  reference  14  are  reproduced 
in  Figure  7.  A  more  in  depth  analysis  is  required  to  predict  spall  size 
and  spall  velocity. 

Kot  et  al_.  ®  studied  spalling  effects  for  large,  thick  concrete 
elements  associated  with  power  plant  structures.  Unlike  the  previous 
described  reference,  the  decay  of  the  peak  pressure  through  the  wall  and 
the  non-uniform  nature  of  the  blast  loading  are  considered.  The  later 
factor  takes  into  consideration  the  added  complication  of  the  spalling 


of  concrete  due  to  normal  and  oblique  incidence,  with  the  determination 

that  no  spall  is  anticipated  for  incidence  angles  greater  than  45°. 

The  set  cf  graphs  available  in  reference  16  is  not  directly  applicable 

to  buried  structures  as  the  blast  load  parameters  for  scaled  distances 
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were  obtained  from  AFM88-22  .  However,  similar  response  would  be  ob¬ 
served  for  buried  structures  with  comparable  loading  functions.  Several 
scaled  graphs  are  presented  which  relate  spall  thickness,  spall  number, 
and  spall  velocity  to  scaled  distance,  wall  thickness,  angle  of  incidence, 
etc.  Some  important  conclusions  can  be  documented  as  a  consequence  of 
the  analysis  of  material  presented  in  reference  16. 

For  walls  where  no  attenuation  of  the  shock  wave  in  the  wall 
is  assumed  to  take  place,  spall  thickness  depends  on  the  assumed  profile 
of  the  compression  blast  wave.  Since  the  blast  wave  is  steepest  at  the 
shock  front  and  becomes  less  steep  as  the  peak  pressure  decays,  the 
first  spall  is  the  thinnest.  However,  it  also  possesses  the  highest 
velocity.  The  authors  also  note  that  ti.e  total  depth  of  spall  is  nearly 
constant,  being  approximately  equal  to  one  half  the  wave  length  of  the 
compression  wave.  Also,  the  angle  of  incidence  has  a  pronounced  effect 
on  spall  thickness  and  *  >  'city.  Kinetic  energy  per  unit  area  for 
oblique  incidences  .  jbstantially  larger  than  for  normal  incidence. 

When  the  wall  thickness  effects  c  z'  decay  of  the  blast  loading  intensity 
is  considered,  the  spall  thickness  increases  with  increasing  wall  thickness 
while  the  spall  velocity  decreases. 

Kot  et  al_.^  also  computed  the  structural  response  in  the  form  of 
wall  displacements  for  impulsive  loads.  A  circular  yield  line  pattern 


is  assumed  to  form  in  a  large  structural  wail.  The  radius  of  this 
yield  pattern  is  determined  by  evaluating  circles  with  incrementally 
larger  radii  until  the  maximum  value  of  deflection  is  determined.  A 
graph  from  reference  16  showing  scaled  concrete  wall  velocities  due  to 
impulsive  air  blast  loading  is  given  in  Figure  8.  This  graph  illus¬ 
trates  the  coupling  between  spalling  and  wall  motion  which  results  in 
high  velocity  ejection  of  most  spall  debris.  It  should  be  pointed  out 
that  modeling  the  structural  response  by  ultimate  strength  yield  line 
procedures  requires  that  the  cross-section  possess  flexural  capacity, 
which  in  reality  is  destroyed  by  the  preceeding  spallation.  This  leaves 
the  member  with  only  the  membrane  capacity  of  the  reinforcing  steel. 

Conclusions 

The  extension  of  the  rigid  plastic  analysis  procedure  to  include  the 
effects  of  axial  constraint  on  the  member's  response  furnishes  an  important 
previously  overlooked  link  necessary  to  explain  a  range  of  experimentally 
observed  failure  modes.  The  extension  strain  consistent  with  axial 
restraint  provides  an  explanation  for  the  experimental  response  observed 
where  the  overall  deflection  is  not  excessive,  but  the  damage  to  the 
concrete  is  extensive.  The  initial  material  degradation  due  to  the 
intensity  of  the  blast  pressure  oscillations  would  destroy  the  flexural 
capacity  of  the  composite  cross-section,  leaving  only  the  membrane 
capacity  of  the  reinforcing  steel  to  determine  the  structural  response. 
While  the  approximate  closed  form  analysis  procedure  presented  demonstrates 
the  significance  of  including  membrane  response,  there  is  definitely  the 
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need  to  develop  a  more  sophisticated  numerical  technique  in  order  to 
examine  the  effects  of  extension  strain  more  closely. 
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APPENDIX  A 


A  review  of  numerous  experimental  test  reports  revealed  significant 
differences  among  test  specimens,  their  restraint  and  loading  conditions,  and 
the  interpretation  of  breach  loading  conditions.  The  later  is  to  be  expected 
considering  the  different  failure  modes  possible  for  reinforced  concrete  slabs 
that  are  subjected  to  blast  pressures  that  range  in  magnitude  from  slightly 
above  the  static  collapse  load  to  several  times  larger  than  f',  the  statically 

V 

determined  cylinder  strength  of  concrete.  Therefore,  no  attempt  is  made  to 
arbitrarily  define  breach  conditions.  Rather,  test  specimens  that  sustained 
severe  damage  are  classified  according  to  their  apparent  mode  of  failure. 

The  failure  test  data  used  in  this  report  is  given  in  Tables  1  and  2. 

The  information  in  Table  1  was  used  in  reference  3  to  obtain  the  previously 
stated  empirical  breach  equation,  while  Table  2  contains  additional  test  data 
presently  available.  From  an  examination  of  the  test  data  and  photographs  of 
these  specimens,  the  following  failure  mechanisms  are  delineated  based  on  the 
type  and  degree  of  damage: 

(1)  Yield  Line,  without  concrete  degradation, 

(2)  Yield  line,  with  concrete  degradation, 

(3)  Membrane,  with  severe  spalling  or  scabbing, 

(4)  Direct  concrete  crushing  and  spalling,  with  membrane  net, 

(5)  Shear,  and 

(6)  Overbreach. 

These  assumed  failure  response  mechanisms  have  the  following  characteristics. 

Yield  Line,  without  concrete  degradation.  This  mechanism  is  characterized  by 
plastic  hinge  formation  at  points  of  maximum  flexural  stress  with  the  sections 
of  the  member  between  hinge  lines  remaining  essentially  intact.  Under¬ 
reinforced  concrete  elements  possess  the  necessary  ductility  for  its  response 
to  be  idealized  by  a  rigid,  perfectly  plastic,  constitutive  relationship. 

This  flexural  failure  is  characterized  by  crushing  of  the  concrete  along  the 
yield  line  in  the  flexural  compression  zone  of  the  crossection.  The  depth  of 
this  concrete  loss  is  not  excessive  because  concrete  slabs  are  usually  con¬ 
siderably  under-reinforced.  With  tied  compressive  reinforcement  present,  the 
member  continues  to  respond  in  this  flexural  mode  until  fracture  occurs  due  to 
excessive  straining  of  the  tension  reinforcement.  For  severe  overloads,  the 
compressive  reinforcement  may  also  fracture,  resulting  in  total  disengagement 


tween  hinges. 


:1eld  Line,  with  concrete  degradation.  In  addition  to  the  previously  described 
flexural  failure  associated  with  the  formation  of  the  plastic  hinge(s), 
concrete  failure  occurs  directly  due  to  the  intensity  of  the  blast  pressure. 

Membrane ,  with  severe  spalling  or  scabbing.  When  the  edges,  or  at  least  the 
corners,  of  a  slab  are  restrained  against  in-plane  motion,  lateral  deflections 
beyond  a  certain  level  (A  >  0.3  t)  are  accompanied  by  stretching  of  the  middle 
surface^0.  As  the  magnitude  of  the  maximum  deflection  increases,  the  resulting 
membrane  force  can  significantly  aid  or  even  predominate  in  carrying  the  lateral 
loads.  The  final  deflected  shape  of  reinforcing  steel  in  test  specimens  failed 
with  this  mechanism  is  consistent  with  the  shape  observed  for  pure  membranes 
subjected  to  distributed  loads.  The  analytical  and  experimental  results  of  ref¬ 
erence  11  support  the  concept  of  this  membtoiie  type  failure  response  for  thin, 
fixed-edge  flat  plates  subjected  to  explosive  pressure  loadings.  Spalling, 
which  is  caused  by  the  shock  pressure  of  a  blast  being  transmitted  through  an 
element,  results  in  a  tension  failure  in  the  concrete  normal  to  its  free 
surface.  Scabbing  is  also  a  tensile  failure  of  the  concrete,  but  is  associated 
with  the  large  strains  In  the  reinforcement  at  the  later  stages  of  ductile 
response  of  a  reinforced  concrete  element.  The  velocities  of  scabbed  fragments 
of  concrete  are  lower  than  the  velocities  of  spalling  fragments. 

Direct  concrete  crushing  and  spalling,  with  membrane.  The  following  description 
from  reference  22  of  the  response  of  a  member  to  a  contact  explosion  also 
describes  the  mode  of  failure  under  consideration.  Near  the  point  of  detonation, 
the  magnitude  of  the  compressive  pulse  is  sufficient  to  crush  the  concrete 
material.  The  pulse  decays  rapidly  as  it  propagates  through  the  member  until 
no  additional  crushing  is  possible  or  the  structure  is  breached.  If  the  later 
does  not  occur,  a  compression  wave  is  reflected  by  the  inner  surface  as  a 
tension  pulse  wave.  If  the  tensile  strength  of  the  concrete  is  exceeded, 
rupture  will  occur  and  the  spall  will  fly  off  in  response  to  the  trapped  mo¬ 
mentum.  Additional  ejection  of  crushed  and  spalled  material  occurs  as  the  tension 
wave  propagates  back  to  its  point  of  detonation.  The  flexural  capacity  of  the 
member  is  lost,  leaving  only  the  membrane  capacity  of  the  reinforcing  net. 


Shear.  This  is  a  catastrophic  failure  mode  in  which  there  is  very  little 
evidence  of  flexural  behavior.  The  affected  area  of  wall  is  typically  blown 
Into  the  structure  nearly  intact,  with  boundaries  showing  direct  shear  failure, 
rather  than  diagonal  tension  failure. 

Overbreach,  Catastrophic  damage  is  obvious,  with  concrete  disentegrated  and 
blown  away.  The  reinforcing  steel  does  not  span  the  opening  and  is  in 
disarray. 

Figure  A1  shows  the  different  slab  support  configurations  used  in  the 
various  testing  programs.  The  experimental  breach  data  is  plotted  in  Figure  A2 
with  each  data  point  symbol  Indicating  the  failure  mechanism  for  that  test 
specimen. 

Except  for  one  test,  the  data  for  X  >  1 .4  consists  entirely  of  the  results 
from  slabs  supported  on  two  legs.  The  slabs  in  these  tests  fall  in  flexural 
through  the  formation  of  yield  lines.  The  author  observed  that  because  of  the 
lack  of  In-plane  translation  restraint  on  the  slab,  the  membrane  capacity  of 
the  slab  is  not  utilized  as  It  might  be  In  a  real  box  structure.  The  distinctive 
yield  line  (plastic  hinge)  Is  not  observed  to  form  In  any  of  the  test  specimens 
which  are  provided  In-plane  restraint  against  translation  by  their  supports. 

Tables  1  and  2  Indicate  that  the  majority  of  the  test  specimens  contain 
considerably  less  steel  reinforcement  than  one  would  expect  to  find  in  actual 
hardened  structures.  No  reference  that  attributes  this  reduced  reinforcement 
value  to  model  scaling  was  found.  While  conservative,  the  use  of  steel  ratios 
significantly  less  than  that  of  actual  structures  should  have  a  pronounced 
effect  on  the  specimens  with  .the  larger  breach  loading  distances  where  structural 
member  response,  rather  than  direct  crushing  of  the  concrete  material,  is  the 
primary  cause  of  failure. 

For  scaled  breach  distances  less  than  1.4  there  appears  to  be  direct 
material  degradation  present  due  to  the  initial  intensity  of  the  blast  pressure 
on  the  restrained  test  specimens.  The  two  test  specimens  in  this  range  that 
did  fail  without  direct  material  crushing  are  unrestrained  two-legged  test 
specimens.  These  two  specimens  also  contain  a  steel  ratio  p  =  0.0025  for  the 
tensile  reinforcing  steel  in  one  direction.  This  definition  for  steel  ratio, 
consistent  with  conventional  reinforced  concrete  practice,  is  equal  to  one- 
fourth  of  the  amount  reported  in  reference  2  as  the  volume  of  steel  in  the 
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specimen.  This  minimal  amount  of  reinforcement,  equivalent  to  what  is  normally 
considered  an  amount  of  steel  simply  for  temperature  and  shrinkage  control, 
probably  accounts  for  the  unexpected  yield  line  failure  mechanism  in  the  high 
intensity  range  of  scaled  breach  distances.  It  should  be  pointed  out  that  the 
yield  line  failure  mechanism  would  appear  essentially  the  same  as  a  modulus  of 
rupture  type  failure  typical  of  a  concrete  flexural  member  that  is  effectively 
unreinforced. 
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Equations  of  Motion  Derived  Using  the  Lagrange  Equation 
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Figure  3.  Strength  Interaction  Diagram  Showing  Both 
Axial  Compression  and  Axial  Tension 


38 


39 


Figure  4.  Simple  beam: (a)  Stationary  deformations  and  forces  (b)  Non-stationary 

deformations  and  forces. 


v  =  Wall  Velocity,  ft/sec 
m  =  Mass  per  Unit  Area,  Slug/ft 
W  =  Charge  Weight,  lb 


I _ I _ I _ I _ 1 _ I _ L 

.2  0.6  1,0  1.4  1.8  2.2  2.6 


Scaled  Standoff  Distance,  ft/lb1^3 


Figure  8.  Scaled  Concrete  Wall  Velocities  Due  to  Impulsive 
Air  Blast  Loading  (From  Reference  16) 


Scaled  truck 


Type  1:  Two  Legged  Support 


Type  2:  Four  Legged  Support 


Figure  Al.  Support  for  Slab  Specimens 


Figure  A2.,  Failure  Mechanisms  for  Blast  Loadings 

43 


wn-Tje. 

r 7.1  mux 


xo 


/o 


/S 


X/ 


4o  \o.oofo 


4o  XAoofo 


4o  I  ooofo 


4.0  I  o.ooso 


4.o 


4.  0  I  oooFo 


4.0  aooSo 


4o  o.ooso 


4.0  \  00050 


A0O63 


0.00  63 


m  <w 


0.03  \>6o\*4r.s 


H  A 

Ot)  M 


o.oo  \*6.o  \*ozr 


333\*6.0  Woo\/.06 


>6.0 \*4X0\aT  I  464 


033  >64  >¥30100  \l14\lS0 


*60  1*4X5100  \#.00\/J9 


133\>6jo\hxs\/.0  \7.14\0S0 


>Lo \»4io\  SO  \J000  \i.st 


*6.0  U IS  to  \x7.oo\t.03 


4.00  \* 6m 


4.oo  \>6o 


>6m 


mm 


JO.ao 

046 

17.40 

A6f 

4.00 

Alt 

LW*  . 
(«/«,*; 

vt- 

0.0*00 

/o.o 

o.oi67 

JO.O 

o.o/ff 

/o.o 

30386 

7.o 

00*33 

to 

O.OXXZ 

i.o 

0.84/3 

4.o 

0. 0037 

/AO 

A  0063 

tl.o 

AOf/6 

6.0 

00370 

60 

O.OI7S 

Af.O 

Table  1.  Experimental  Failure  Data  for  Breach  Equation  of  Reference  3. 
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